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Abstract. — We describe a new software package that may be used to determine stellar and atomic parameters 
by matching observed spectra with synthetic spectra generated from parameterized atmospheres. A nonlinear least 
squares algorithm is used to solve for any subset of allowed parameters, which include atomic data (log gf and van der 
Waals damping constants), model atmosphere specifications (T e ff, log 3 ), elemental abundances, and radial, turbulent, 
and rotational velocities. LTE synthesis software handles discontiguous spectral intervals and complex atomic blends. 
As a demonstration, we fit 26 Fe I lines in the NSO Solar Atlas (Kurucz et al. 1984), determining various solar and 
atomic parameters. 

Key words: radiative transfer — methods: numerical — techniques: spectroscopic — stars: fundamental 

parameters — sun: fundamental parameters — atomic data 


1. Introduction 

Traditional stellar spectroscopy may be divided into four 
steps: spectroscopic observation, data reduction, measure- 
ment of stellar parameters (e.g. rotational velocities, ef- 
fective temperature, abundances), and astrophysical inter- 
pretation. Large telescopes, high resolution spectrographs, 
and digital detectors have made it possible to routinely 
obtain spectra at a resolution of 200 000 with noise below 
0.3% for stars down to 8th magnitude. Reduction software 
has kept pace with improvements in instrumentation, but 
we have yet to take full advantage of modern comput- 
ers in the measurement phase of the analysis. There is 
no reason to discard the rich information content in mod- 
ern spectra, simply because traditional analyses only use 
equivalent widths. By matching synthetic spectra to the 
observed spectrum, we can minimize errors due to blends 
and continuum placement. In addition, modern computers 
make it possible to search parameter space for the model 
which best matches the observations. This approach al- 
lows one to characterize errors in each parameter and de- 
cide whether it is better to fix a parameter to an externally 
determined value or to solve for the parameter using the 
observed spectrum as a constraint. 

In this paper, we describe a new software package 
(Spectroscopy Made Easy, SME) which allows one to find 
the best fit to an observed spectrum for a given set of 
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assumptions. Remaining discrepancies may then be inter- 
preted as inaccuracies in atomic line parameters and/or 
inadequacies in the adopted model atmosphere (e.g. gran- 
ulation is ignored). SME can be used to determine atomic 
parameters (log gf and van der Waals damping) for indi- 
vidual transitions and/or global stellar parameters (Kad^ 
Vsin i, log#, abundances, V^cro, and Knacro)- SME 
may be divided conceptually into three main components: 
a spectrum synthesis library written in C++, parame- 
ter optimization code written in IDL and intended for 
batch execution, and a widget-based user interface which 
is also written in IDL. Analysis tools with related goals, 
but somewhat different approaches are being developed by 
Hill (Hill et al. 1996) and Cowley (Cowley 1996). Below, 
we describe how SME works and illustrate its use with an 
analysis of the solar spectrum. 


2. Radiative transfer 

Spectral synthesis is the most computationally intensive 
part of SME, so we use a very fast algorithm adapted from 
the program SYNTH (Piskunov 1992). Currently in SME, 
we assume: local thermodynamical equilibrium (LTE), 
plane parallel geometry, no bulk flows, no molecular line 
opacity, and negligible magnetic field. SME handles depth 
dependent microturbulence that is (unlike SYNTH) the 
same for all lines. 
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The spectral synthesis code is implemented as a C+ + 
library, which is dynamically linked with IDL at execution 
time. This allows us to combine the speed of C++ with 
the high level functionality and user interface of IDL. The 
C++ library contains modules for computing ionization 
balance, continuous opacities, emergent intensities, and 
contribution functions, as well as I/O modules for com- 
municating with IDL. 

2 . 1 . Data input 

Before the C++ library can calculate an intensity spec- 
trum, four types of information must be specified: a model 
atmosphere, elemental abundances, atomic line data, and 
a set of wavelength intervals. Because of the perturbative 
nature of the nonlinear least squares algorithm described 
below, SME can often make use of certain results from 
preceding calculations. Accordingly, each type of data re- 
quired by the C++ library can be modified independently 
from IDL, allowing a partial respecification of a radiative 
transfer request. The radiative transfer code checks that 
all required data have been supplied and then calculates 
any quantities that have been invalidated by new request 
specifications. For example, changing the atomic param- 
eters for one line has no effect on ionization equilibrium 
and affects the intensity spectra in only one wavelength 
interval. Adding a line from a previously unrepresented 
element, however, will force a recalculation of the ioniza- 
tion balance and all subsequent steps. 


Continuous absorption coefficients (in cm 2 g _1 ) are 
computed at the endpoints of each wavelength interval, 
using the ionization equilibrium depth grid. Continuous 
opacity routines were adapted from the ATLAS 9 spec- 
trum synthesis code (Kurucz 1993). 


2.3. Solution of radiative transfer equation 


Specific intensities are determined by solving the integral 
form of the radiative transfer equation: 


/a 



e- T ^S A (r)dr, 


( 2 ) 


where r is the optical depth and S\(r) is the source func- 
tion (i.e. the Planck function, B\ y given our assumption of 
LTE). Equation (2) is solved using Gaussian quadrature: 


h = (Ti)-Wi, (3) 

1=1, N 


where r t and w t are the quadrature nodes and weights (see 
Table 1) for a basis set of Hermite polynomials. We find 
that 10th order quadrature gives fractional integration er- 
rors less than 2 10“ 5 for realistic model atmospheres. 


Table 1. Nodes and weights of 10th order quadrature formula 
used for solving the radiative transfer equation 


Ti 


Wi 


2.2. Ionization equilibrium and continuous opacity 

Ionization equilibrium is determined by solving the Saha 
equation: 


N, __ (2ttA:T) 3 /2 2Ui f AE\ 

N t - 1 “ h 3 N e ' Ui - 1 <>XP V kTJ' 


( 1 ) 


which gives the ratio of two consecutive ionization states 
as a function of electron density and temperature. For the 
partition functions, £/, we use polynomial approximations 
which describe the first six ionization stages of the first 99 
elements in the periodic table (Kurucz 1979). 

Ionization fractions are computed on a finer depth grid 
than that of the original model atmosphere because ioniza- 
tion can change more rapidly than other quantities. The 
number of subdivisions is determined at each depth to 
ensure that the the fractional change of each ionization 
state is less than 10%. Normally, only a few depth inter- 
vals need to be divided into a maximum of 20 subintervals 
each in order to satisfy this condition. The temperature 
and electron density in each subinterval is determined by 
parabolic interpolation in the original grid. Ionization bal- 
ance is determined for hydrogen, helium, and every other 
element present in the line list. Neutral hydrogen and he- 
lium populations are used later for computing eollisional 
line broadening. 


0.137793470540 
0.729454549503 
1 808342901740 
3.401433697855 
5.552496140064 
8.330152746764 
11.843785837900 
16.279257831378 
21.996585811981 
29.920097012274 


3.084411 15765e-01 
4.011 19929155e-01 
2. 18068287612e-01 
6.20874560987e-02 
9.501 51697518e-03 
7.53008388588e-04 
2.82592334960e-05 
4.24931398496e-07 
1.83956482398e-09 
9.91182721961e-13 


The model atmosphere is tabulated either as a function 
of optical depth (r re f) at a reference wavelength (typically 
at 5000 A) or mass column density (pdx). For each wave- 
length point in the spectrum, B\ must be calculated at the 
optical depths, r 7 . used in the quadrature. Conversion be- 
tween optical depth, r\, and the model depth parameter is 
accomplished by solving one of two first order differential 
equations: 


<1 Vr 

A'ref 

(4) 

11 

I<\ + Ax ’ 

pdx 

1 

(5) 

(l dn 

A a + A a 


where either Eqs. (4) or (5) is used, depending on whether 
T re f or pdx is the independent depth parameter. In these 
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equations, K\ and I< re f are the continuous opacities (in 
cm 2 g _1 ) at A and A re f, and A\ is the line opacity given 

by 



A 2 

y/ncAXp 



H(a,v) 


( 6 ) 

where g , /, and N are the statistical weight (degeneracy) 
of the lower level, the oscillator strength, and the number 
density of the atoms in the lower level. The summation is 
over all lines contributing at given wavelength. The Voigt 


function, 


H(a , v) 



exp (— 3/ 2 ) 
(v — y) 2 + a 2 


d y 


( 7 ) 


is evaluated as a function of depth, using the Humlicek 
approximation (Humlicek 1982). 

A 6th order Runge-Kutta technique with appropriate 
boundary conditions is used to integrate Eqs. (4) or (5) 
and find the value of the independent depth parameter 
at each of the quadrature nodes. Integration is done on 
the depth grid constructed in the ionization equilibrium 
calculation (Sect. 2.2). This approach has two major ad- 
vantages over direct integration of Eq. (2). First, r re f is 
a smooth, monotonic function of r* that is well sampled 
on the finer depth grid. Second, the integration progresses 
downwards from the surface, so it can be truncated when 
r\ > 1 . 

For each Runge-Kutta step, we interpolate the main 
atmospheric variables: temperature, number densities of 
electrons and atoms, ionization fractions for contributing 
lines, and microturbulence. We calculate Voigt profiles, 
including the effects of radiative damping, the quadratic 
Stark effect, and van der Waals broadening. Damping con- 
stants (r rad , r 4 , r 6 ) that are not provided with the line 
list are approximated using formulas from Unsold (1955), 
Cowley (1971), or Griem (1968). Stark broadening in hy- 
drogen and helium must be treated as a special case. For 
the B aimer series of Hydrogen, we use the approximation 
of Vidal et al. (1971), and for 20 optical lines of He I, 
we use tables of Stark broadening and central wavelength 
shifts calculated by Dimitrijevic (1971). 

Synthetic spectra are computed on an adaptive wave- 
length grid to produce uniform precision when interpo- 
lating the spectrum. We calculate the specific intensity 
at the endpoints of each wavelength interval and at each 
line center. We then insert new wavelength points midway 
between existing points until the accuracy of the inter- 
polated profile is better than a selectable threshold (see 
Vincent et al. 1993 for more details). This algorithm re- 
sults in more wavelength points where the spectrum has 
higher curvature. 

The ultimate precision of the synthetic spectrum at an 
arbitrary wavelength is controlled by two factors: the pre- 
cision of the radiative transfer solution at each point in 


the irregularly spaced wavelength grid, and the interpola- 
tion accuracy between these grid points. We can achieve a 
fractional precision of 10“ 4 , which is more than adequate 
in most situations, given current uncertainties in model 
atmospheres and line formation mechanisms. 


3. Fitting the observations 

3.1. Atomic parameters 

SME requires line data for all atomic transitions of inter- 
est. Mandatory data are element name, ionization state, 
wavelength, excitation energy of the initial state, and 
log (//. Default approximations are used to estimate the 
radiative, Stark, and van der Waals damping constants, 
if they are not explicitly specified. SME initially reads 
atomic data from one or more text files in the format re- 
turned by a “stellar” request to the Vienna Atomic Line 
Database (VALD, Piskunov et al. 1995) electronic mail 
server. VALD critically selects atomic line data for those 
transitions that contribute significantly to the spectral 
synthesis. Once loaded into SME, atomic data may be 
saved in an internal format for later use. 

SME can be used to solve for log gf values and/or van 
der Waals damping constants. Parameters may be varied 
individually on a line by line basis or globally via a cor- 
rection to all van der Waals damping constants. Global 
enhancement of the van der Waals damping predicted by 
Unsold (1955) is common practice, for example when mod- 
eling neutral iron lines (Holweger et al. 1991). Using SME, 
atomic parameters may be adjusted to best match an ob- 
served spectrum of some star (the Sun, for example). The 
new “astrophysical” atomic parameters are often an im- 
provement over theoretical values or low quality labora- 
tory data, and can be used in subsequent analyses of other 
stars. 

3.2. Model parameters 

In order to synthesize a spectrum, SME requires a model 
atmosphere. The required atmospheric data are tempera- 
ture, electron number density, atomic number density, and 
mass density. These data are specified on a depth grid pa- 
rameterized either by integrated mass column or by con- 
tinuum optical depth at a reference wavelength (5000 A, 
for example). SME could easily be modified to calculate 
(rather than read) densities, as long as radiation pressure 
is negligible. 

Model atmospheres are often tabulated on a grid of 
{T ef r, log </, [Fe/H]}. SME uses bicubic spline interpola- 
tion of the logarithm of each atmospheric variable to gen- 
erate an atmosphere for any desired {T’ e flr> log p, [Fe/H]} 
contained within the limits of the grid. The interpolated 
model is constructed on the mass column (or continuum 
optical depth) scale of the closest model in the grid. Grid 
interpolation makes it easier to generate an atmospheric 
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model for a specific {T eff , log#, [Fe/H] } or to solve for the 
{r eff , log#, [Fe/H]} that give the best agreement with an 
observed spectrum. 

3.3. Disk integration of specific intensity 

Specific intensities must be integrated over the stellar sur- 
face to produce a flux spectrum. Techniques for doing 
this vary in both speed and accuracy. For nonrotating 
stars, specific intensity depends only on p, making one- 
dimensional Gaussian quadrature accurate and efficient. 
For rotating stars, annular regions of constant p span a 
range of rotational velocities, and a two dimensional sur- 
face integral is formally required to compute the flux. To 
increase the speed of this calculation, one may first cal- 
culate the flux profile ignoring rotation and then convolve 
the flux profile with a broadening kernel designed to mimic 
the effects of rotation. This approximation is inadequate 
for precise spectroscopy of stars rotating at less than about 
15 km s -1 (Gray 1992). 



LoSiot^sooo) 

Fig. 1 . Comparison of the temperature distributions for the 
HM, SME, and Kurucz solar models. Arrows indicate the 
depths of formation for the cores of each of the 26 Fe I lines in 
our analysis 

In SME we employ a recently developed disk integra- 
tion technique in which the star is divided into annu- 
lar regions, each characterized by a single intensity spec- 
trum. The effects of solid-body rotation are handled by 
convolving each intensity spectrum with an analytic ker- 
nel that properly weights the velocity shifts in the corre- 
sponding annulus. In this manner, calculation of the flux 
spectrum is reduced to a one-dimensional sum of con- 
volutions, resulting in good speed and excellent accuracy 
(Valenti & Anderson 1995). By default SME uses seven p 
angles, which gives an integration accuracy better than 
0.1% in most circumstances. The number of p angles can 
be changed to vary the speed and accuracy of calcula- 
tion. The effects of isotropic Gaussian inacroturbulenee 


are treated by a single convolution of the rotationally 
broadened flux profile. When radial and tangential macro- 
turbulent. velocities are unequal, each annulus in the disk 
integration has its own macroturbulent broadening kernel. 

3.4. Matching the observations 

SME simplifies the determination of atomic and stellar 
parameters by using an observed spectrum as a model 
constraint. When empirically matching an observed spec- 
trum, the behavior of the spectrograph and detector must 
be included in the model. An instrumental profile is speci- 
fied in tabular form or as an analytic function (a Gaussian, 
for example). The synthetic spectrum is spline interpo- 
lated onto a uniform wavelength grid with points spaced 
finely enough to preserve all information contained in both 
the model spectrum and the instrumental profile. The in- 
strumental profile is interpolated onto (or generated on) 
the same wavelength grid and convolved with the syn- 
thetic spectrum. 

Each segment of synthetic spectrum is independently 
shifted in wavelength and linearly scaled in amplitude un- 
til x 2 is minimized. For each trial wavelength shift, scaling 
and slope adjustment of the model spectrum are found di- 
rectly by linear least squares matching of observed “con- 
tinuum” points. If necessary, the synthetic spectrum is 
integrated over each observed pixel to mimic wavelength 
binning by the detector. Continuum points are specified 
explicitly by the user or determined automatically by 
SME, based on how close the points lie to the maximum 
flux in the spectral segment. The presence of spectral lines 
in the “continuum” fitting regions only affects continuum 
placement if the depth of model lines is wrong. 

To find the optimum wavelength offset, we calculate 
X 2 for a sequence of wavelength shifts (and correspond- 
ing continuum adjustments). A parabola is fit to \ 2 as a 
function of wavelength shift. The parameters at the min- 
imum of the parabola are then used to match the syn- 
thetic spectrum to the observations. The wavelength ad- 
justments measure and correct errors in radial velocity, 
transition wavelengths of the dominant lines, and zero 
point errors in the observed wavelength scale. The con- 
tinuum adjustments mimic wavelength variations in spec- 
trograph throughput and correct for errors in observed 
continuum placement and/or model continuous opacities. 

3.5. Parameter optimization 

As indicated in Sects. 3.1 and 3.2, SME can solve for 
various atomic and/or model parameters, using an ob- 
served spectrum as a constraint. Free parameters are opti- 
mized simultaneously, using a \ 2 minimization algorithm 
described by Marquardt (Marquardt, 1963; Press et al. 
1986). Sequential optimization of different sets of parame- 
ters may be done by running SME once for each parameter 
set. 
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Fig. 2. A comparison of 25 Fe I observed and model line profiles after adjusting solar parameters and individual oscillator 
strengths The top of each panel shows the synthetic spectrum (solid line) and the observations (filled circles for points fit by 
SME- open circles for points that were ignored). The bottom shows residuals multiplied by 10, where the thicker line corresponds 
to the solid circles above. Horizontal dashed lines indicate ±1 % error levels. Panels are sorted by (and labelled with) excitation 
potential. Central wavelengths (in A) are also given. Wavelength offsets (in A) and residual intensity for all lines are indicated 

in the bottom left hand panel 
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The Marquardt algorithm mixes a traditional gradient 
search with approximate parabolic minimization of the 
X 2 surface near the minimum. Approximate second par- 
tial derivatives in the curvature matrix are constructed 
from first partial derivatives, which in turn are found by 
recomputing synthetic spectra at small offsets in each free 
parameter. For each set of partial derivatives, three to five 
model calculations are required to choose a good mix be- 
tween a gradient search and parabolic minimization. The 
new parameters are used to recalculate partial derivatives 
and the process repeats until \ 2 changes by less than some 
fractional amount, typically 0.2%. 

Convergence to a fractional limit of 0.2%, typically re- 
quires only two to five major iterations (for which new 
partial derivatives are recalculated). As a rule of thumb, 
we find that model convergence requires about 3.5 N -f 13 
model spectrum calculations. The inverse of the final cur- 
vature matrix is the covariance matrix, which contains for- 
mal errors for each parameter and degeneracies between 
each pair of parameters. 

4. Solar example 

To illustrate the capabilities of SME, we solved 
for certain solar and atomic parameters by match- 
ing synthetic spectra to a flux spectrum of the Sun 
(Kurucz et al. 1984). We concentrated on a set of 26 
Fe I lines with precisely measured laboratory oscilla- 
tor strengths (Blackwell et al. 1995, BLS and references 
therein). The equivalent widths of these lines were used 
by BLS to determine the abundance of iron in the Sun. 
For demonstration purposes, we decided to repeat their 
analysis, using SME to match hue profiles. The results 
are shown in Figs. 2 and 3. After completing our tests, 
we learned that Holweger et al. (1995, HKB) strongly dis- 
agree with various interpretations in BLS. It is not our 
goal here to try and resolve this debate, but we do make 
a few specific comments below, based on our tests with 
SME. 

4.1. Observational and atomic data 





Fig. 3. The same comparison as in Fig. 2 for two of the lines. 
The residuals are shown for laboratory log gf (dashed line) and 
astrophysical log gf determined by SME (solid line) 


We used a continuum normalized spectrum of the entire 

solar surface (Kurucz et al. 1984) as input data for SME. tervals with expected fractional depths greater than 0.2% 
The noise in this FTS spectrum is 0.1 to 0.2%, with a spec- bi the Sun. For the 26 Fe I lines used in the abundance 
tral resolving power of 5.2 10 5 . We ignored instrumental analysis, we adopted Oxford \og(gf) values as tabulated 
broadening, which is negligible since the instrumental pro- m BLS, except that log (gf) = -5.040 (not -5.004) for the 
file is one-tenth as wide as the spectral lines. Blends in the line at 6574.24 A. Damping constants from Kurucz CD- 
wings of the 26 Fe I lines were identified by comparing an ROM 18 wen 1 initially provided by VALD, though the van 
initial model spectrum with observations. Blends which <ler Waals damping constants were subsequently increased 
were poorly modeled due to inaccurate or missing atomic by SME, as described below, 
data were marked and ignored in subsequent analyses. 

Atomic data were obtained from VALD for 5 A wide A n n * . . . 

spectral intervals containing the 26 Fe I lines studied by 

BLS. We used VALD “stellar” requests to obtain the best SME makes it relatively easy to solve for any combination 
available atomic, data for o39 lines in these spectral in- of stellar and atomic parameters, but it doesn^ simplify 
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the decision of which parameters to vary and which to hold 
constant. The constraints provided by the observed spec- 
trum lead to a certain accuracy for each derived param- 
eter. If a particular parameter is known more accurately 
from previous work, then the parameter should be set to 
the known value. Since accuracy (unlike precision) can be 
difficult to assess, the decision about which parameters to 
fix is often subjective. 


Table 2. Derived parameters of the solar models 


Teff 

(K) 

Fe/H 

1 micro 

(kin/s) 

Fm aero 
(km/s) 

V sin t 
(km/s) 

af 6 

RMS 

(%) 

5734 

-4.50 

0.71 

2.12 

1.55 

2.75 

0.95 

5731 

-4.51 

0.74 

2.11 

1.63 

(2.5) 

0.96 

5775 

-4.46 

0.79 

2.18 

1.50 

(2.0) 

1.00 

5808 

-4.42 

0.83 

2.19 

1.55 

(1.5) 

1.05 


For our solar example, we set the surface gravity to 
log g— 4.44 cm s -2 and used VALD to fix log gf and the 
relative values of Tq. Initially, we let AITe vary with our 
other free parameters after which we fixed ATe to partic- 
ular values of interest. We then simultaneously varied the 
five remaining parameters: T e ff, log(Fe/H), F m icro, Knacro, 
and V sin i. Although the actual solar T e ff is well deter- 
mined, we let this parameter vary in our models in order to 
empirically determine the theoretical Kurucz model that 
best reproduces the observed spectrum. Rotation period 
as a function of solar latitude is also known, but the best 
choice for V sin i in a uniformly rotating model is unclear, 
so we solve for this parameter as well. 

Table 2 gives solar parameters that we derived using 
SME. Iron abundances in column two are given as loga- 
rithmic ratios with respect to the total number of atoms. 
Parentheses around entries in column six indicate that 
A T 6 was set to the tabulated value. The last column gives 
RMS residuals for points deeper than 1% of the line depth. 
Formal uncertainties are all much less than the tabulated 
precision, but systematic errors and correlated uncertain- 
ties are certainly larger. 

Model convergence for the A T 6 = 2.5 case required 
68 spectrum synthesis operations, including calculation 
of eight separate partial derivatives with respect to five 
free parameters. Reduced \ 2 improved from 5800 to 740, 
changing by only 0.05 in the final iteration. The large value 
of x 2 implies that no combination of model parameters 
adequately reproduces a high quality solar spectrum. The 
line RMS of 0.96% is partially due to errors in log gf, but 
the line shapes are also incorrect. In particular, observed 
line profiles are asymmetric (probably due to granulation), 
whereas our static, one-component models produce sym- 
metric model profiles. The continuum RMS of 0.14% is 


well above the noise, due to line asymmetries and uniden- 
tified weak blends. 

Table 3. Astrophysical log gf values determined by SME with 
associated atomic data and line-by-line fit quality 


A 

(A) 

Xex 

(eV) 

Lab 

logs/ 

Solar 

logs/ 

Transition 

Line 

RMS 

4389.245 

0.052 

-4.583 

-4.562 

a 5 D 3 - z 7 FS 

1.26 

4445.471 

0.087 

-5.441 

-5.446 

a 5 D 2 - z 7 F! 

1.07 

5247.049 

0.087 

-4.946 

-4.936 

a 5 D 2 - z 7 DS 

0.76 

5250.208 

0.121 

-4.938 

-4.886 

a 5 Do - z 7 D? 

0.76 

5701.545 

2.559 

-2.216 

-2.149 

b 3 F2 4 - y 3 D| 

0.54 

5956.692 

0.859 

-4.605 

-4.589 

a 5 F 5 - z 7 P 4 

0.99 

6082.708 

2.223 

-3.573 

-3.614 

a 5 Pi - z 3 P? 

0.48 

6136.993 

2.198 

-2.950 

-2.995 

a 5 P 2 - y 5 D? 

0.79 

6151.617 

2.176 

-3.299 

-3.333 

a 5 P 3 - y 5 DS 

0.50 

6173.341 

2.223 

-2.880 

-2.898 

a 5 Pi - y 5 Dg 

0.63 

6200.314 

2.609 

-2.437 

-2.393 

b 3 F2 2 - y 3 F° 3 

0.77 

6219.279 

2.198 

-2.433 

-2.492 

a 5 P 2 - y 5 D| 

0.61 

6265.131 

2.176 

-2.550 

-2.583 

a 5 P 3 - y 5 Bl 

0.46 

6280.616 

0.859 

-4.387 

-4.402 

a 5 F 5 - z 7 F? 

1.08 

6297.792 

2.223 

-2.740 

-2.774 

a 5 P i -y 5 D? 

0.65 

6322.690 

2.588 

-2.426 

-2.384 

b 3 F2 3 - y 3 F° t 

0.69 

6481.869 

2.279 

-2.984 

-2.962 

a 3 P2 2 - y 5 D? 

0.60 

6498.937 

0.958 

-4.699 

-4.677 

a 5 F 3 - z 7 F° 3 

0.71 

6574.225 

0.990 

-5.004 

-5.030 

a 5 F 2 - z 7 FS 

0.30 

6593.871 

2.433 

-2.422 

-2.396 

a 3 H 5 - z 5 Gg 

0.94 

6609.110 

2.559 

-2.692 

-2.664 

b 3 F2 4 - z 3 G 4 

0.60 

6625.021 

1.011 

-5.336 

-5.381 

a 5 Fi - z 7 F? 

0.21 

6750.150 

2.424 

-2.621 

-2.594 

a 3 P2 i - z 3 P? 

0.55 

6945.203 

2.424 

-2.482 

-2.480 

a 3 P2 i - z 3 P| 

0.47 

6978.850 

2.484 

-2.500 

-2.474 

a 3 P2 0 - z 3 P? 

0.45 

7723.205 

2.279 

-3.617 

-3.521 

a 3 P2 2 - z 3 D! 

0.50 


4.3. Comparison with expected values 

Before comparing individual parameters in Table 2 with 
expected values from the literature, we should emphasize 
again that parameters have correlated uncertainties and 
are subject to systematic errors. Each parameter has a 
unique effect on line shape, but no combination of pa- 
rameters exactly matches the observed profiles. While at- 
tempting to match the observed profiles, parameters de- 
viate from their true values by an unknown amount to 
partially compensate for deficiencies in the model. Such 
errors are difficult to quantify. 

One way to explore possible systematic errors is to try 
a range of reasonable values for one parameter and solve 
for the remaining free parameters. The results of such an 
experiment are shown in Table 2, where we set A Tq to 2.5, 
2.0, and 1.5, covering the range of values deemed accept- 
able by Holweger et al. (1991). The remaining five param- 
eters were all varied to achieve the optimum fit. The fit 




